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Abstract 

In this paper, we calculate neutron and proton decay rates and branch- 
ing ratios in a predictive SO(IO) SUSY GUT which agrees well with low 
energy data. We show that the nucleon lifetimes are consistent with the 
experimental bounds. The nucleon decay rates are calculated using all 
one-loop chargino and gluino dressed diagrams regardless of their chiral 
structure. We show that the four-fermion operator Cjk{uRdjj^){dkLi^TL), 
commonly neglected in previous nucleon decay calculations, not only con- 
tributes significantly to nucleon decay, but, for many values of the ini- 
tial GUT parameters and for large tan /3, actually dominates the decay 
rate. As a consequence, we find that rp/r„ is often substantially larger 
than the prediction obtained in small tan/3 models. We also find that 
gluino-dressed diagrams, often neglected in nucleon decay calculations, 
contribute significantly to nucleon decay. In addition we find that the 
branching ratios obtained from this realistic SO (10) SUSY GUT differ 
significantly from the predictions obtained from "generic" SU(5) SUSY 
GUTS. Thus nucleon decay branching ratios, when observed, can be used 
to test theories of fermion masses. 



1 Introduction 



It can easily be seen that Grand Unified Theories [GUTS] contain baryon num- 
ber violating operators that produce proton decay Q. However, determining 
the decay rate is a much more involved task. In supersymmetric [SUSY] GUTS, 
the nucleon decay amplitude is directly proportional to the inverse of an effective 
color triplet Higgs mass Mt (resulting from effective dimension 5 baryon number 
violating operators) ^ |^ multiplied by a product of the Yukawa couplings of 
this color triplet Higgs to quarks and leptons. Thus, obtaining a theoretical 
prediction for nucleon decay rates depends critically on two factors: 

1. obtaining bounds on the effective color triplet mass, and 

In this regard, it has been noted that Mt also affects the prediction for 
as through threshold corrections at Mgut- Thus, bounds on Mt can be 
obtained via the experimental constraint on |@j ^ but only if one 
has a complete SUSY GUT, valid above Mqut- 

2. predictions for the relevant Yukawa couplings. 

The nucleon decay branching ratios are sensitive to these color triplet 
Higgs-quark-quark and Higgs-quark-lepton Yukawa couplings. These cou- 
plings are completely determined in any predictive theory of fermion masses 
and mixing angles, but they are typically NOT identical to the Yukawa 
couplings responsible for quark and lepton masses. 

In addition, one must take these effective dimension 5 operators and renor- 
malize them from Mgut to Mz- Then, at the weak scale, effective dimension 
6 operators are obtained by closing the squark and/or slepton lines into a loop 
via chargino or gluino exchanges. Hence the decay rate depends sensitively on 
soft SUSY breaking parameters. Finally the effective dimension 6 operators 
are renormalized from Mz to the nucleon mass and then an effective chiral 
Lagrangian analysis is used to obtain lifetimes and branching ratios. 

In a recent paper ^ [paper I], we proposed several complete SO(IO) SUSY 
GUTs, valid above Mgut- At Mgut, these models had the desired feature that 
they reproduced the effective fermion mass operators of model 4 of Anderson, 
et al. [01 [paper II]. Moreover, we found that if we required that our models 
contain in their superpotentials all terms not forbidden by the symmetries of the 
respective models, one and only one additional effective fermion mass operator 
was generated for two of the models, while no additional effective fermion mass 
generating operators were generated in a third model. We showed that for 
reasonable values of the GUT scale parameters, we can obtain values of as 
consistent with experiment, with the effective color triplet Higgs mass that 
enters into proton decay rates Mt as large as 10^^ GeV or even larger. We 
argue, however, that the scale ^ 10^^ GeV is, in fact, a natural upper bound 
for Mt. 
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Note, the actual color triplets' masses are of order the GUT scale. Such a 
large value for Mt in comparison to the GUT scale is obtained in our model 
using the Dimopoulos-Wilczck mechanism for doublet-triplet splitting In 
Appendix 1 , we review how the Dimopoulos-Wilczek mechanism can be used to 
obtain Mt > Mqut- 

In this paper, we present the details of our calculations of the nucleon decay 
rates for one of our models, referred to as model 4(c), which includes an addi- 
tional "13" mass operator. In a forthcoming paper, Blazek, Carena, Wagner, 
and one of us (S.R.) will show that model 4(c) fits all low energy data to within 
la, while model 4 (without the additional operator) agrees with all data at the 
2cr level only §. Our main results for the predictions of the proton and neutron 
decay rates and branching ratios in model 4(c) are found in Tables AlA, AlB 
and A6. 

In addition we have studied the sensitivity of our predictions to different 
factors. We have compared the predictions for our model 4(c) to those of models 
4(a) through (f) of Anderson et al. This tests the sensitivity of the predicted 
branching ratios to the quality of the fit for ferniion masses and mixing angles. 
We find that branching ratios can differ by factors as large as 10. 

We have also compared the predictions of our models with those of small 
tan/3 minimal SUSY SU(5) GUTs. We find that some results in the large tan/3 
regime are qualitatively different than for small tan f3. For example, certain 
dimension 6 operators with the chiral structurc|^ LLRR tend to dominate over 
their LLLL counterparts for fJ-{Mz)/mi/2 > 1. In this limit, the ratio Tp/rn, for 
example, is sensitive to the quality of the fcrmion mass fit and is substantially 
larger than what it is in small tan/3 SUSY GUTs. 

Finally we discuss the sensitivity of our predictions to neglecting either the 
gluino or chargino exchange diagrams. We find both contributions to be signif- 
icant. 

The paper is organized as follows — in sections 2, 3 and 4 we discuss the 
calculational ingredients. Namely, in section 2, we discuss GUT models and the 
physics from Mqut to Mz] in section 3, SUSY loops at Mz and the resulting 
dimension 6 baryon violating operators; and in section 4, the physics from Mz to 
the nucleon mass and a summary of the numerical procedure and our results. In 
section 5, we discuss our results and the aforementioned sensitivities to different 
factors. 



^ LLLL, LLRR, and RRRR refer to four-fermion operators pairing four left-handed Weyl 
fermions, pairing two left handed Weyl fermions and the conjugates of two right-handed Weyl 
fermions, and pairing the conjugates of four right-handed Weyl fermions, respectively. Specific 
examples of each type of four-fermion operator can be found in Table 3, infra. 
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2 The low energy effective operators generating 
nucleon decay, and their renormaUzation 

Recall that in paper II, it was shown that four effective fermion mass operators, 
denoted ©33, ©23, C12, and O22, could be used to fit fermion masses and mixing 
angles reasonably well. For model 4 of that paper, the first three of these 
operators were 

©33 = I63 lOi I63 

O23 = I62 ^ lOi ^ I63 
A A 

0,2 = I61 (A)3 (_i)3 

<->M 

where I61, I62, and I63 are 16 representations containing the first, second, and 
third generations of fermions of the Standard Model, respectively; lOi is a 10 
representation containing the two Higgs doublets of the Minimal Supcrsymmet- 
ric extension of the Standard Model; Ai, A2, and A arc 45s which get vcvs in 
the B-L, hyperchargo, and X (SU(5) invariant) directions, respectively; and Sm 
is an SO(IO) singlet getting a vev of order Mpianck- There were six different 
choices for the operator O22, labeled a through f. 

022 = 



(a) 


162^ 

<->M 


Ai 
101^162 
A 


ib) 


162=1 
A 


lOi 4^ I62 


(c) 


A 

I62 ^ 

<->M 


Ai 

lOi ^ I62 

<->M 


id) 


I62 lOi 


Ai 
A 


(e) 


I62 lOi 


2 I62 

OM 


(/) 


I62 lOi 


AiSg 
A^ 



where Sg is an SO(IO) singlet getting a vev of order Mqut- In addition, from 
paper I, we discovered that when we built a complete GUT, valid up to energies 
around Mpianck and we allowed all terms in the superspace potential consis- 
tent with the symmetries of the theory, an additional fermion mass generating 
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operator, O13, was generated for certain versionsg of model 4. 



(a) I61 (— )3 lOi (-4) I63 

(b) 

(c) 161 (A)3 (A) 

At low energies, these five operators produce effective operators which gen- 
erate the fermion masses and which are responsible for baryon-number violating 
nucleon decay. 

033 + O23 + O22 + O12 + Ol3 

HuQYjj + HdQYdD + HdLYeE + Q^CggQT + QcqiLT + Uc^dDT + Uc^eET 

HuQYjJ + HdQYdD + HaLY^E + -^Q^Cq^Q QcqiL + j^UcudD UcueE 

with T and T being the color triplet Higgses from lOi. Yu-, Yd, Ye, Cqq, Cqi, Cudi 
and Cue are flavor matrices which can be expressed in terms of seven independent 
real parameters. At the GUT scale, the values of these matrices are 



Yu^ 



Yd = 



Y, 



' 

c 


c 



-IB 






-27C 


-27C 

Ee^'f" 

-IB 


UdDe'^ \ 

I 




-27C 


-27C 

iEe^'t' 
2B 






C 


C 

VqqEe'^ 

\B 


UqqDe'^ 

\B 
A 



-27C UqiDe 

cqi = I -27C yqiEe'^ \B 
^'qiDe^' \B A 



■5 - 



^Note, models d, e and f have the second family I62 coupled directly to lOi and a heavy 
16. If this coupling is as large the third generation Yukawa coupling, then we would obtain 
excessively large flavor changing neutral current processes, such as ^ e + 7. 
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Table 1: y Clebsches for each version of model 4 



model 


Vqq 


Vql 


Uud 


Hue 


a 


-3/4 


3/4 


-5/4 


-3/4 


b 


-3/2 


5/2 


1/2 


-3/2 


c 


-1/2 


3/2 


-1/2 


-1/2 


d 


3/2 


3/2 


-1/2 


3/2 


e 


1/2 


5/2 


1/2 


1/2 


f 


9/4 


3/4 


-5/4 


9/4 



Table 2: u Clebsches for models 4(a) through (c) 



model 


Uu 


< 


Ud 




Ue 


< 


Uqq 


Uql 




Uud 


Kd 






a 


-4/3 


1/3 


-2 


-9 


-54 


3 


1/3 


3 


-9 


-2 


36 


2 


-4/3 


b 









































c 


-4/3 


1/3 


2/3 


-9 


-54 


-1 


1/3 


-1 


-9 


2/3 


36 


2 


-4/3 





-27C 

c„d = I -27C VudEe"^ 
IB 

C 

VueEe'^ -4B I (1) 
-2B 

with the values of the Clebsches given in Tables |l| and |. 

These matrices need to be renormalized from the GUT scale to the elec- 
troweak scale. Due to the no-renormalization theorems of supcrsymmetry, only 
wavefunction renormalizations enter into the calculation of the renormalization 
group equations of these matrices. The RGEs for the matrices are 

- -^[(.YuYUY,Yj)cqq + Cqq{Y^Y^+Y,Yjf (2) 



dt 167r2 

(ysl + 3.g2 + Y^dDcqq] 

dCql _ 1 



[{YuY^+Y,Yj)cqi+CqiiY,Y,^)^ 



dt 167r2 
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dt 167r2 
with t = log(/i/Mz). 

3 Nucleon decay formulas 

At low energies, the matrices Y^, Yd, and Ye are diagonahzed by unitary matrices 
called S and T defined so that 



SuYuTu — Ac — Yi 




and so forth, with all diagonal entries in Yd, and Yf. being real and positive. 
The weak eigenstate basis for the fermionic fields will be represented by primed 
fields, and the mass eigenstate basis will be represented by unprimed letters. 
Additionally, we choose an unprimed basis for neutrinos such that 

S>' = 

so that the lepton-lepton-weak boson vertices are flavor diagonal in the unprimed 
basis.0 

Squarks and sleptons mass matrices are diagonahzed using 6 by 6 matrices 
r defined 



f2, VL = u,d,e 



where $7 is a six dimensional vector of mass eigenstates and fl' and il' are weak 
eigenstates. Cf. notation of Additionally, we define Tq ^ and ^ to be 
6 by 3 matrices so that Tq ^ consists of the first three columns of Fj^ and Fq 
consists of the last three columns. In block matrix notation, 

= ( ^n,L I ^n,R ) 

Since there is no left-handed right-handed neutrino mixing, we will define 
Fj^ so that 



•^If left-handed neutrinos have a mass, the v^, v^, and Vt thus defined are generally not 
mass eigenstates. Since we sum decay rates over all neutrino species in our analysis, the 
convention for defining the neutrino basis is irrelevant for our results. In paper I, the 3 left- 
handed neutrinos remain massless, but this result is easily changed with the introduction of 
GUT scale majorana neutrino masses for heavy singlet neutrinos. 
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with v being a mass eigenstate basis for the sneutrinos. A right-handed 
matrix will not be defined. 

Finally, the chargino mass matrices are diagonalized by matrices U+ and C/_ 
defined 



and 




where Xi 2 mass eigenstates. 

In this notation, various Feynman vertices relevant to nucleon decay are 
given in Appendix 2. 



3.1 Gluino diagrams 

Figure shows the two diagrams which contribute to the four-fermion operator 
Cljl^i'"^ ^^\ufdj){d1vi)eaf3'y Calculating these two diagrams, we find 

^{ud){di^)[G] _ 4 1 „2{t^ F* F F* ^ i'lj';^ t(~ ~ ~i \ 

^ijkl — 3 leTT^Mt^Sl U,L U,L \i' ^ D,L pj^ D,L pj'^qq '^ql "'■g-' li/j "A, «pj 

^D,L p.'i^kL pf^DX ak^DX .k'^,,'^'' c,i''^'rngl{~g, dp, d.)} 

where a, /?, and 7 are color indices, i, i', j, j' , k, and / are fermion flavor indices, 
p and A are squark flavor indices, c^q — SuC^gSj , c^i — SdC^iSj with and 
being Cqq and Cqi, respectively, renormalized to Mz, and 

^ ^ mglogmg ^ mllogml ^ mglogmg 

{ml - ml){ml - m^J {ml - ml) {ml - ml) {ml - ml) {ml - ml) 

where ma equals the mass of particle a, etc.0 

By similar calculations, we can construct the remaining four-fermion opera- 
tors listed in Table ^. Note, we also use Table ^ to define our notation for these 
four-fermion operators. The Cijkis for these operators are 

^{ud){ue)[G] _ 4 1 2fp F* F T* r ^' r M'- m - J ( h 1~I ^ d \~ 

'-'ijkl — Q,„ 93,\^ U,L U,L D,L pj^ D.L p]'(^qq '^ql ^g^yg^Ux^ap) 

^U,L Xi^UX Xi'^UX ak^UX ak'^qq^ ^ Cql^ ^'mgl{g, Ux, Ua)} 

^(dd)iu^)lG] _ 4 1 2p F* r r* r '[j'r J d) 

'^ijkl — 2 -|^g^2^ ^3^^ D,_R, pjJ^ Z5,L pi'J^ Z5,i? (Tfe^ _D,L crfc'Cgg Cq/ mgl [g , Up, ) 

*We use a notation for antisymmctrization on tensor indices in which A^j/.^^^ is equal to 

A-ijklmn ~ ^iZfcjmni ^T^d SO forth. 
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Table 3: Table of all gluino-dressed four fermion operators relevant to nucleon 
decay 



operator type 


LLLL 


LLRR 


RRRR 










^ijkl K )Ke;)ea/37 





^iud){d,y)[G] _ 4 1 2p r* r r* r "-^^ r ^^^'-m-Tlh il^ d\ 

'^iikl ~Qir. -) ir,r U,R\i^ U,L D,R pj^ D.L pj''^qq ^ql ">-g^\yT U-X^apj 

ry(ud)(ue)[G\ _ ^ 1 2p T* T T* r [''^'r .''^'■m-Kn -Ti ^ d ) 

^ijkl — ^ 93'^ U,B. Xi'^ U,L Xi''^ D,R pj'^ D,L pj'^^qq '^ql ^g^\9,Ux,apj 

^iud)(ue)[G] _ 4 1 Op p* p p* [j'j'-;* k]l T(~ ~ 7n 

"-"ij-fei — o,„ ^ 9?,^ U,L Xi^ U.BXi'^ D,L pj^ D,R pi'^^ud ^ue "*gJ 15: "A, Cpj 

•J IVTT iVIt 

^(ud){ue)[G] _ 4 1 2 fp p* p p* -* [I'j'-;* k]l T/~ -- J \_ 

-r(7,fl Ai-rt/,fl Ai'T^^^ o-fc-Tt/.fl o-fc''^Jld'* ■'^ue*^ '''71g/(5, Ua, Uct)} 

where c,; = S^c^iS^ , Cud = T^c^d^d, and Cne = T^c^e^e with cf^ and c^^ 
being c„d and c„e, respectively, renormalized to M^. 

As observed in ref. the contribution of gluino dressed operators is 

zero in the limit that all squarks are degenerate at the electroweak scale. In 
this analysis, however, squarks and sleptons are assumed degenerate at Mqut 
and as a consequence of renormalization group running they are explicitly non- 
degenerate at the weak scale. Thus we retain the gluino contribution in our 
analysis. 

In principle, we could also have four-fermion operators like C'l"^!^''''*'^'"'^' x 

{u*°'u*f,''){dj ei)eap-f pairing two up type quarks in a Weyl index contracted pair. 
However, since we are only interested in the first generation of up quarks, the 
portions of these operators that would be relevant to nucleon decay are identi- 
cally zero. 
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3.2 Chargino diagrams 

Using the vertices in Appendix 2, the chargino diagrams can be readily com- 
puted. According to analyses by many other authors on proton decay in SUSY 
GUTs, the dominant operators would be expected to be C^j^i^'^'^^^^^ {ufdj){dli'i)eai3j 
and Clj^i^^'^^^^\ufd^){u1ei)ea0'y The Cijkis for these operators are equal to 

X [{92Tu,lU+ In - ^U,R^u U+ 2n)VKM]\i 
'>^^u,L\i'^*D,Lpj' cj Cql ''^''m^^ I{Xn ,Ux,dp) 
+ [{92^U,lU+ In - ^U,R^u U+ 2n)VKM]\k 

^[{92^E,lU- In - ^E.rYc U- 2n)]pl 

X [{92T^U,lU+ In - ^U,rYu U+ 2n)VKM]\j 
><'^UX \i''^DX pj'Cqq^''^' ^ql''^^^Xr^IiXn,Ux,dp) 

+[{92Td,lU- m - Td^^Y^ 2n)vlj^U 

X92T, XlU+ inT*D,L pk'K Xl'Cq,^^c/^^'m^J{Xn, dp, vx)} 

The Cijkis for the LLRR and RRRR operators are 

cgT'^^"^^ = ^^^{K 2n{^D,LVUYu UUl 2nirU,LVKMY, )x, 

X^U,L Xi'^D,L pj'cj"^' c/'^^m^J{xn,ux,dp) 

+ [{92^U,lU+ In - ^U,rYu U+ 2n)VKM]xk 
X [i92TE,LU- In - ^E,Rye U- 2n)]pl 
X^U,R Xk'^*E,R pl'Kid^^ KiJ' ''T^Xr^HXn, Ux, ep)| 
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^[{9-2^U,lU+ In - ^U,rYu U+ 2n)VKM]\j 
U,R D,i?. pj'^ud ^ue '^Xrz^KXn, UX, "pj 

4 Numerical procedure and results 

In calculating the following nucleon decay rates, we make the standard universal- 
ity assumptions about the soft SUSY parameters at Mqjjt, except that we allow 
non-universal values for Mh^ and Mh^ -0 We define Mqut such that ai = a2 = 
otGUT at Mgut and define £3, representing the contribution of GUT scale thresh- 
old corrections to the gauge coupHngs, to be (aslMauT) — aQUT)/oiGUT- The 
dimensionless (dimensionful) parameters are renormalized at two (one) loops to 
Mz using the rcnormalization group equations of Martin and Vaughn [p^ , ex- 
cept that Cqq, Cqi, Cud, and Cue are renormalized at one loop using the equations 
of this paper. 

Rcnormalization of the Cijkis from Mz to 1 GeV is taken into account 
by multiplying them by the Al calculated in and then chiral Lagrangian 

^Note that if the messenger scale of SUSY breaking is Mpianck then our analysis is not 
completely self-consistent. In any complete SUSY GUT defined up to an effective cut-off scale 
M > Mq, the interactions above Mq will renormalize the soft breaking parameters. This 
will, in general, split the degeneracy of squark and slepton masses at Mq even if they are 
degenerate at M. On the other hand, bounds on flavor changing neutral current processes, 
severely constrain the magnitude of possible splitting. Thus these corrections must be small. 
In addition, in theories where SUSY breaking is mediated by gauge exchanges with a messenger 
scale below (but near) Mq, the present analysis is expected to apply unchanged. Since in 
this case squarks and sleptons will be nearly degenerate at the messenger scale. The Higgs 
masses, on the other hand, are probably dominated by new interactions which also generate a 
H term. It is thus plausible to expect the Higgs masses to be split and independent of squark 
and slepton masses. The parameter Aq could also be universal at the messenger scale. 
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techniques IQ are used to obtain nucleon decay amplitudes. Formulas for 
nucleon decay rates in terms of the chiral Lagrangian parameters are contained 
in Appendix 3. 

The decay rates depend heavily on the chiral Lagrangian factors a and /3 
where 

(3Uik) = e„0^<O|(u"d^)w'^|proton(k)>, 

Q;C/(k) = eQ/37<0|(u*"c?*'^)u'''|proton(k)> 

and U{\i) is the left handed component of the proton's wavefunction. See, e.g., 
ref. jl|L It is known that |/3| = \a\ ^ ^ and that |/3| ranges from .003 to .03 
GeV*^ |17|. Lattice calculations have not reduced the uncertainty in lattice 
calculations have reported \(3\ as low as .006 GeV'^ and as high as .03 GeV'^ 
[|8|. Additionally, the phase between a and (3 is not widely reported, although 
a relatively recent lattice calculation suggests that P = —a [0. Therefore, we 
have left the phase between a and f3 a free variable, and report three values for 
many of the quantities predicted in our tables. Namely, the max (min) referred 
to in the tables is the value for the quantity predicted when the phase between 
a and /3 is such that the quantity is maximized (minimized). Hence, each entry 
in the max and min columns uses a different value of arg(/3/a). 

In the following tables, we have calculated decay rates using Mt — 10^^ GeV 
and = .003 GeV^. In paper I, we showed that without any fine tuning this 
value of Mt can be made consistent with the measured value for as. In addition, 
we argued that it seems unnatural to have Mt much bigger than 10^^ GeV in 
the sense that in order to have Mt much bigger than 10^^ GeV, there would 
need to be a supermassive electroweak doublet in the GUT desert with mass 
many orders of magnitude lower than the GUT scale itself and at least an order 
of magnitude lighter than any other particle getting mass around the GUT 
scale. Therefore, the values presented in the tables are roughly upper bounds 
on the nucleons' lifetimes, based on naturalness^ Since all decay rates scale as 
^ 003^ GeV^ 1/'^°'^ -^^' ^^'^ lifetimes for different values of Mt and \P\ can easily 
be extracted from the tables. 

In Tables Al through A6 [the A tables], we calculate nucleon decay using 
model 4(c) (including the O13 operator of paper I), since it appears to be the 
model which best fits the low energy data . Initial values for the dimensionless 
Yukawa parameters; soft SUSY parameters; and Mqut, olgut, and tan/3 are 
taken from the global analysis of Blazek et al. This global analysis 
shows that these values of the GUT parameters are consistent with electroweak 

note of caution - it was also shown in ref. that chiral Lagrangian techniques 

overestimate the amphtude eQ,^^<7r''|(n°'d'^)n^|p> by at least a factor of 2.4 (for (3 = 0.006). 
Thus the proton decay rate is overestimated by almost a factor of 6 or more in this case. As a 
result, for any given |/3|, the actual nucleon lifetimes for that /3, could be a factor of 6 or more 
larger than the results reported here, as extrapolated from our tables for that value of |/3|. 
These remarks are indicative of the theoretical uncertainty in the calculation due to strong 
interaction effects. 
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symmetry breaking and the experimental bounds on the sparticle masses; and 
that, for the particular values of the soft SUSY breaking parameters used, these 
parameters give the best global fit to 20 low energy observables, including ex- 
perimental measurements for the gauge couplings; fermion masses and mixing 
angles; and b sj, with per degree of freedom < Tables AlA and AlB 
contain the lifetimes for p K^V and n — > K^v and compares these rates with 
the rates of decay into the other significant decay modes involving spin zero 
mesons, for various values of the GUT scale parameters. These are the main 
results of this paper. Note however that only the three or four most significant 
decay modes are included in these tables. 

In the subsequent tables we evaluate the relative contributions to these rates 
from different sources - (LLRR vs. LLLL operators) or (gluinos vs. charginos). 
Tables A2A, A2B, ASA, and A3B compare the contributions of chargino and 
gluino diagrams to the total decay rate. Table A4 compares the contribution of 
LLLL versus LLRR operators for decays into anti-neutrinos, for each of the three 
anti-neutrino species. Tables ASA and A5B compare the relative importance of 
each generation of anti-neutrino to the total proton decay rate, for decay modes 
involving anti-neutrinos. Finally, Table A6 contains the values of the GUT scale 
parameters used in Tables Al through A5. 

Tables Bl through B6 [the B tables] contain information similar to that in 
the A tables except that they compare models 4(a) through (f), without the 
Oi3 operator. We used values for the initial (GUT scale) parameters, taken 
from unreported data from the collaboration of ref . |^ , which are consistent 
with electroweak symmetry breaking and the experimental bounds on sparticle 
masses, and which give predictions agreeing to within 2.1(7 with experimental 
measurements for gauge couplings, fermion masses and mixing angles, and h 
S"f. This comparison gives us information on the model dependence of nucleon 
decay branching ratios. Note that models 4(a) through (f) (without the ©13 
operator) give identical results for fermion masses and mixing angles. This is 
because the contribution of the different 022 operators to the 22 entry of the 
Yukawa matrices all give the same Clebsch relation 0:1:3 for u:d:e matrices 0. 
They however have different Clebsch relations for the 22 entry of the matrices 
Cqq, Cqi, Cud, and Cue relevant for nucleon decay. 

In addition, the comparison of model 4(c), with the O13 operator, to models 
4(a) through (f), without the O13 operator, gives us information on the sensi- 
tivity of predictions for nucleon decay with respect to the quality of the fit for 
fermion masses and mixing angles. In order to compare the runs in the A tables 
directly with the runs in the B tables, we chose the GUT scale parameters such 
that for each run in the B tables, the values for a^^y, Mqut, £3; the soft SUSY 
breaking parameters tan/3, ^, toi/2, toq: '^Hdi and Aq; and the Yukawa 

parameter A are nearly the same as they are for the run of the same name in 
the A tables. 
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run 
no. 


t{p - 


* K+V)/(10 


32yrs) 


r(p^7r+I7) 

r(p-»if+i7) 


r(p^ifV+) 
r(p^/f+i7) 


X 10^ 


r(p^7rV+) 

r(p->K+T7) 


X 10^ 


r(p 
r(p 




X 10^ 




max 


13= -a 


min 


max 


13 = -a 


mm 


max 


/^=- 


a mm 


max 


f3 = —a min 


max 


(3 = —a min 


I 


26. 


14. 


14. 


0.46 


0.39 


0.39 


0.52 


0.29 


0.29 


0.28 


0.16 


0.16 


0.099 


0.055 


0.055 


II 


60. 


37. 


37. 


0.37 


0.33 


0.32 


0.53 


0.33 


0.32 


0.29 


0.18 


0.18 


0.10 


0.062 


0.061 


iii(i) 


220. 


130. 


98. 


1.1 


0.74 


0.69 


0.31 


0.19 


0.14 


0.12 


0.073 


0.055 


0.028 


0.017 


0.013 


111(2) 


150. 


97. 


74. 


1.5 


1.1 


0.94 


0.37 


0.24 


0.18 


0.11 


0.071 


0.054 


0.017 


0.011 


0.0084 


111(3) 


110. 


76. 


58. 


1.5 


1.1 


1.0 


0.34 


0.23 


0.18 


0.092 


0.063 


0.049 


0.011 


0.0078 


0.0060 



Table AlA: Partial mean lifetime for proton decaying into kaon plus anti-neutrino and 
ratios of the rates of proton decay into various decay products versus rate of decay 
into kaon plus anti-neutrino for various values of the GUT scale parameters, when the 
Ci3 operator is included. 



run 

no. 


T(n ^7^°F)/(103Vs) 


r(n^7r"l7) .^,2 


r(n^)7i/) ■,f.2 


r(n^7r .|f^2 


max P = —a min 


max P = —a min 


max P = —a min 


max P = —a min 


I 

II 
iii(i) 

111(2) 
111(3) 


3.8 2.3 2.3 
12. 7.3 6.9 
14. 12. 9.8 

6.4 5.8 5.0 

4.5 4.2 3.7 


3.3 3.2 3.2 
3.3 3.2 3.2 
3.5 3.5 3.4 
3.2 3.1 3.1 
3.2 3.1 3.1 


0.37 0.17 0.17 
0.49 0.25 0.23 
0.13 0.088 0.035 
0.079 0.055 0.022 
0.063 0.045 0.019 


0.083 0.051 0.050 
0.11 0.070 0.067 
0.016 0.014 0.011 
0.0094 0.0085 0.0073 
0.0076 0.0070 0.0061 



Table AlB: Partial mean lifetime for neutron decaying into kaon plus anti-neutrino 

and ratios of the rates of proton decay into various decay products versus rate of decay 
into kaon plus anti-neutrino for various values of the GUT scale parameters, when the 
Oi3 operator is included. 



run 




pchargino 






pCHargmo 

%^^+- 






pCtiargino 




no. 




ptotal 






ptotal 






ptotal 






max 


P — —a 


min 


max 


P = —a 


min 


max 


P — ~a 


min 


I 


1.1 


1.1 


0.90 


1.1 


1.1 


0.87 


0.99 


0.99 


0.99 


II 


1.3 


1.2 


0.82 


1.3 


1.2 


0.73 


1.0 


1.0 


1.0 


iii(i) 


1.9 


1.4 


0.45 


1.4 


1.3 


0.71 


1.2 


1.2 


1.2 


111(2) 


1.7 


1.3 


0.52 


1.3 


1.2 


0.79 


1.1 


1.1 


1.1 


111(3) 


1.6 


1.3 


0.57 


1.2 


1.2 


0.83 


1.1 


1.1 


1.1 



Table A2A: Ratios of the rate of proton decay that would occur if chargino diagrams 

contributed only versus total proton decay rate for the three most dominant decay 
modes, for various values of the GUT scale parameters, when the ds operator is 
included. 
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run 
no. 



I 
II 

iii(i) 

111(2) 
111(3) 



otal 



max (3 = —a 



mm 



1.1 
1.3 
1.4 
1.2 
1.2 



1.1 0.87 

1.2 0.73 

1.3 0.71 
1.2 0.80 
1.1 0.84 



rto 



max f3 = —a 



mm 



1.1 
1.3 
1.4 
1.3 
1.2 



1.1 0.87 

1.2 0.73 

1.3 0.71 
1.2 0.79 
1.2 0.83 



chargii 



max /3 = —a 



mm 



1.1 0.99 0.97 

1.2 0.99 0.96 
2.2 0.44 0.40 
2.2 0.42 0.39 
2.2 0.43 0.41 



0.81 0.81 0.81 

0.64 0.64 0.64 

0.43 0.43 0.43 

0.55 0.55 0.55 

0.63 0.63 0.63 



Table A2B: Ratios of the rate of neutron decay that would occur if chargino diagrams 
contributed only versus total neutron decay rate, for various values of the GUT scale 
parameters, when the ds operator is included. 



run 




pgluino 






pglu.n^_ 






pgluino 




no. 




ptotal 






ptotal 






ptotal 






max 


(3= -a 


min 


max 


13 = -a 


min 


max 


j3 = —a 


min 


I 


0.019 


0.010 


0.010 


0.0060 


0.0038 


0.0038 


0.00038 


0.00038 


0.00038 


II 


0.085 


0.053 


0.051 


0.036 


0.024 


0.023 


0.00067 


0.00067 


0.00067 


iii(i) 


0.25 


0.15 


0.11 


0.035 


0.030 


0.023 


0.020 


0.020 


0.020 


111(2) 


0.16 


0.10 


0.076 


0.016 


0.014 


0.012 


0.013 


0.013 


0.013 


111(3) 


0.10 


0.071 


0.055 


0.0099 


0.0090 


0.0078 


0.010 


0.010 


0.010 



Table ASA: Ratios of the rate of proton decay that would occur if gluino diagrams 
contributed only versus total protou decay rate for the three most dominant decay 
modes, for various values of the GUT scale parameters, when the Oiz operator is 
included. 



run 




■pgluino 






pgluino 






pgluino 






pgluino 




no. 




ptotal 






ptotal 






ptotal 






ptotal 






max 


[3 = —a 


min 


max 


(3 = —a 


min 


max 


(3 = —a 


min 


max 


13 = —a 


min 


I 


0.0055 


0.0034 


0.0034 


0.0060 


0.0038 


0.0038 


0.025 


0.025 


0.019 


0.0099 


0.0099 


0.0099 


II 


0.033 


0.021 


0.020 


0.036 


0.024 


0.023 


0.11 


0.11 


0.087 


0.043 


0.043 


0.043 


Iii(i) 


0.033 


0.028 


0.023 


0.035 


0.030 


0.023 


0.83 


0.41 


0.32 


0.20 


0.20 


0.20 


111(2) 


0.014 


0.012 


0.010 


0.016 


0.014 


0.012 


0.60 


0.28 


0.22 


0.20 


0.20 


0.20 


111(3) 


0.0087 


0.0080 


0.0070 


0.0099 


0.0090 


0.0078 


0.47 


0.22 


0.17 


0.16 


0.16 


0.16 



Table A3B: Ratios of the rate of neutron decay that would occur if gluino diagrams 
contributed only versus total neutron decay rate for various values of the GUT scale 
parameters, when the C13 operator is included. 
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run 

no. 


/pLLRR 
/ p^/f + -, 


/ pLLRR 


\ / pLLLL 

V r>^K+77, 




77, 77^1 77 j- 




I 

II 
iii(i) 

111(2) 
111(3) 


0.000084 0.011 2.0 
0.000086 0.011 1.8 
0.00021 0.029 3.8 
0.00026 0.040 4.6 
0.00031 0.048 5.7 


0.00019 0.024 7.3 
0.00020 0.025 6.4 
0.00056 0.068 10. 
0.00067 0.097 14. 
0.00080 0.12 17. 



Table A4: Ratios of the rate of proton decay that would occur if LLLL operators 
contributed only versus the rate of proton decay that would occur if LLRR operators 

contribxitcd only, for each of the three anti-neutrino generations, for various values of 
the GUT scale parameters, when the O13 operator is included. 



run 




LLLL 




LLRR 


LLLL 

jf+-p 


LLRR 

Tp ^ 


LLLL 

Tp ^ 


no. 




LLRR 


\ 


LLRR 


LLRR 


LLRR 


LLRR 

Tp ^ 


I 


0.082 


6. 


i X 10-6 


1.5 


0.017 


0.49 


II 


0.10 


8.6 X 10-6 


1.9 


0.021 


0.56 


iii(i) 


0.035 


7.2 X 10-6 


0.60 


0.017 


0.26 


111(2) 


0.032 


8.2 X 10-6 


0.49 


0.020 


0.22 


111(3) 


0.026 


8.1 X 10-6 


0.41 


0.020 


0.18 



Table A5A: Ratios of partial decay rates for p — > K'^V, which compare the importance 
of the LLLL and LLRR operators for each generation of anti-neutrino versus contri- 
bution of the LLRR operator of the third generation anti-neutrino for various values 
of the GUT scale parameters, when the O13 operator is included. 



run 


LLLL 

Tp ^ 




LLRR 

Tp ^ 

,r + -e 




LLLL 

Tp ^ 
-+-p 


LLRR 

Tp ^ 


LLLL 

Tp 


no. 


LLRR 


^ 


LLRR 
Tp ^ 


^ 


LLRR 


LLllk 


LLRR 


I 


0.026 




\ X 10- 


6 


0.48 


0.012 


0.14 


II 


0.030 


6.1 X 10- 


6 


0.59 


0.015 


0.16 


iii(i) 


0.0095 


5.3 X 10- 


6 


0.18 


0.013 


0.098 


111(2) 


0.0081 


5.4 X 10- 


6 


0.13 


0.013 


0.071 


111(3) 


0.0067 


5.4 X 10- 


6 


0.11 


0.013 


0.058 



Table A5B: Ratios of partial decay rates for p — » tt 17, which compare the importance of 
the LLLL and LLRR operators for each generation of anti-neutrino versus contribution 
of the LLRR operator of the third generation anti-neutrino for various values of the 
GUT scale parameters, when the da operator is included. 
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run 
no. 


I 


II 


III(l) 


111(2) 


111(3) 




24.43 


24.36 


24.51 


24.65 


24.75 


Mqut 


2.498 X 10^*^ 


3.172 X 10^^ 


3.327 X 10^*^ 


2.857 X 10^*5 


2.513 X lO^'' 


£3 


-0.04760 


-0.04886 


-0.04342 


-0.04420 


-0.04550 


A 


0.7640 


0.8067 


0.8523 


0.8867 


0.8872 


B 


0.05259 


0.05439 


0.05630 


0.05882 


0.05956 


C 


0.0001096 


0.0001155 


0.0001213 


0.0001231 


0.0001226 


E 


0.01251 


0.01308 


0.01360 


0.01397 


0.01397 




1.066 


1.041 


1.020 


1.023 


1.038 


D 


0.0004633 


0.0004944 


0.0005064 


0.0005691 


0.0005665 


5 


5.698 


5.706 


5.698 


5.742 


5.744 


tan j3 


52.77 


54.38 


55.39 


55.86 


55.92 


fi{Mz) 


80.0 


80.0 


160. 


240. 


300. 


mi/2 


280. 


240. 


170. 


170. 


170. 


Too 


400. 


700. 


1400. 


1400. 


1400. 




706.4 


994.6 


1858. 


1859. 


1855. 




635.6 


865.3 


1599. 


1591. 


1585. 


Ao 


322.2 


458.4 


-982.4 


-1079. 


-1274. 



Table A6: Values of the GUT scale parameters used in Tables Al through A5. All 
dimensions in GeV units. 
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run 


o 


T 




r(p^7r+I7) 


r(p^ifV+) 


r(p^,rV+) 




no. 


d 
1 




1032yrs 


T{p^K+V) 


r{p^K+v) 


T{p^K+u) 


r{p^K+v) 




a 




7.3 


0.22 


0.0020 


0.0011 


0.00039 




b 




3.5 


0.17 


0.0032 


0.0017 


0.00060 


I 


c 




11. 


0.25 


0.0012 


0.00067 


0.00024 




d 




4.1 


0.15 


0.0047 


0.0023 


0.00082 




e 




14. 


0.25 


0.0028 


0.0013 


0.00045 




f 




2.0 


0.13 


0.0049 


0.0024 


0.00086 




a 




23. 


0.21 


0.0028 


0.0015 


0.00054 




b 




10. 


0.16 


0.0040 


0.0021 


0.00076 


II 


c 




36. 


0.24 


0.0018 


0.00098 


0.00034 




d 




10. 


0.15 


0.0050 


0.0025 


0.00087 




e 




39. 


0.23 


0.0036 


0.0016 


0.00055 




f 




4.9 


0.13 


0.0051 


0.0026 


0.00090 




a 




59. 


0.30 


0.00078 


0.00034 


0.00011 




b 




68. 


0.28 


0.0023 


0.0012 


0.00039 


iii(i) 


c 




68. 


0.31 


0.00052 


0.00017 


0.000044 




d 




32. 


0.24 


0.0016 


0.00072 


0.00025 




e 




57. 


0.30 


0.00080 


0.00025 


0.000076 




f 




21. 


0.21 


0.0021 


0.00097 


0.00034 




a 




26. 


0.31 


0.00048 


0.00017 


0.000048 




b 




33. 


0.32 


0.0013 


0.00059 


0.00019 


111(2) 


c 




29. 


0.32 


0.00039 


0.000094 


0.000019 




d 




18. 


0.27 


0.0010 


0.00042 


0.00014 




e 




26. 


0.32 


0.00053 


0.00014 


0.000036 




I 




13. 


0.24 


O.OOll 


0.()()0(i2 


0.00021 




a 




18. 


0.31 


0.00040 


0.00012 


0.000031 




b 




23. 


0.34 


0.0010 


0.00040 


0.00012 


111(3) 


c 




20. 


0.32 


0.00034 


0.000073 


0.000012 




d 




14. 


0.28 


0.00081 


0.00031 


0.000099 




e 




19. 


0.32 


0.00044 


0.00010 


0.000024 




J 




11. 


().2() 


0.0011 


0.00047 


O.OOOIG 



Table BIA: Partial mean lifetime for proton decaying into kaon plus anti-neutrino and 
ratios of the rates of proton decay into various decay products versus rate of decay 
into kaon plus anti-neutrino for various values of the GUT scale parameters, when the 
Oi3 operator is not included. For all entries, j3 = —a. 
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run 


o 




r(n^7r"l7) 


r(n— *r/i/) 


r(n^7r-/i+) 


no. 


d 
1 


1032yrs 


T(n^K°V) 


T(n^KOV) 






a 


2.4 


0.036 


0.0020 


0.00071 




b 


1.4 


0.033 


0.0047 


0.0013 


I 


c 


3.1 


0.036 


0.00095 


0.00038 




d 


1.9 


0.035 


0.010 


0.0022 




e 


4.0 


0.036 


0.0039 


0.00073 




f 


1.0 


0.035 


0.012 


0.0025 




a 


8.0 


0.035 


0.0028 


0.0011 




b 


4.2 


0.033 


0.0063 


0.0018 


II 


c 


11. 


0.036 


0.0015 


0.00059 




d 


4.8 


0.035 


0.010 


0.0024 




e 


12. 


0.036 


0.0044 


0.0010 




f 


2.5 


0.034 


0.011 


0.0026 




a 


16. 


0.039 


0.0011 


0.00018 




b 


17. 


0.036 


0.0038 


0.00060 


iii(i) 


c 


17. 


0.038 


0.00075 


0.000082 




d 


9.9 


0.037 


0.0015 


0.00044 




e 


14. 


0.036 


0.00037 


0.00012 




f 


7.4 


0.036 


0.0026 


0.00068 




a 


6.6 


0.039 


0.00066 


0.000086 




b 


7.4 


0.037 


0.0020 


0.00027 


111(2) 


c 


7.0 


0.038 


0.00051 


0.000045 




d 


5.0 


0.037 


0.00085 


0.00023 






6.1 


0.036 


0.00028 


0.000064 




f 


4.1 


().();-!7 


0.0015 


().()()():-!<s 




a 


4.5 


0.039 


0.00054 


0.000060 




b 


5.1 


0.037 


0.0014 


0.00018 


111(3) 


c 


4.8 


0.038 


0.00044 


0.000034 




d 


3.7 


0.037 


0.00064 


0.00017 






4.3 


0.037 


0.00026 


0.000047 




f 


3.2 


O.U3cS 


U.OOll 


U. 00028 



Table BIB: Partial mean lifetime for neutron decaying into kaon plus anti-neutrino 
and ratios of the rates of neutron decay into various decay products versus rate of 
decay into kaon plus anti-neutrino for various values of the GUT scale parameters, 
when the da operator is not included. For all entries, /3 = —a. 
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run 


m 




pcliargino 
p—*K+~ 






■pcliargino 






pcliargino 
p^K^ fj.+ 




no. 


d 
1 




ptotal 






ptotal 






ptotal 

p^K^ fj,'^ 








max 


P = -a 


min 


max 


/3 = -a 


min 


max 


(3 = -a 


min 




a 


1.2 


1.1 


0.85 


1.1 


1.1 


0.87 


1.0 


1.0 


1.0 




b 


1.2 


1.2 


0.93 


1.2 


1.2 


0.86 


1.0 


1.0 


1.0 


I 


c 


1.2 


1.1 


0.83 


1.1 


1.1 


0.89 


0.99 


0.99 


0.99 




d 


1.2 


0.99 


0.98 


1.2 


0.92 


0.91 


1.0 


1.0 


1.0 




e 


1.2 


0.91 


0.91 


1.1 


0.91 


0.91 


1.0 


1.0 


1.0 




f 


1.2 


1.0 


1.0 


1.2 


0.97 


0.94 


1.0 


1.0 


1.0 




a 


1.5 


1.3 


0.72 


1.4 


1.3 


0.72 


1.0 


1.0 


1.0 




b 


1.4 


1.4 


0.90 


1.4 


1.4 


0.75 


1.0 


1.0 


1.0 


II 


c 


1.4 


1.3 


0.65 


1.3 


1.2 


0.75 


1.0 


1.0 


1.0 




d 


1.4 


1.0 


1.0 


1.4 


0.88 


0.86 


1.0 


1.0 


1.0 




e 


1.4 


0.84 


0.84 


1.3 


0.81 


0.81 


1.0 


1.0 


1.0 




f 


1.4 


1.1 


1.1 


1.4 


0.97 


0.93 


1.0 


1.0 


1.0 




a 


1.9 


1.5 


0.53 


1.5 


1.3 


0.69 


1.1 


1.1 


1.1 




b 


2.6 


2.4 


0.34 


2.0 


2.0 


0.47 


1.2 


1.2 


1.2 


iii(i) 


c 


1.6 


1.4 


0.62 


1.3 


1.2 


0.76 


1.2 


1.2 


1.2 




d 


2.7 


0.48 


0.42 


2.0 


0.53 


0.50 


1.1 


1.1 


1.1 




e 


1.7 


0.65 


0.65 


1.3 


0.76 


0.76 


0.98 


0.98 


0.98 




f 


2.8 


0.53 


0.42 


2.3 


0.50 


0.42 


1.1 


1.1 


1.1 




a 


1.5 


1.3 


0.65 


1.3 


1.2 


0.78 


1.1 


1.1 


1.1 




b 


2.2 


2.0 


0.42 


1.6 


1.6 


0.59 


1.2 


1.2 


1.2 


111(2) 


c 


1.4 


1.2 


0.72 


1.2 


1.2 


0.83 


1.1 


1.1 


1.1 




d 


2.2 


0.52 


0.47 


1.6 


0.62 


0.60 


1.0 


1.0 


1.0 




e 


1.4 


0.74 


0.74 


1.2 


0.83 


0.83 


0.99 


0.99 


0.99 




f 


2.5 


0.51 


0.41 


1.9 


0.56 


0.50 


1.1 


1.1 


1.1 




a 


1.4 


1.3 


0.70 


1.2 


1.2 


0.81 


1.1 


1.1 


1.1 




b 


1.9 


1.8 


0.48 


1.5 


1.5 


0.65 


1.1 


1.1 


1.1 


111(3) 


c 


1.3 


1.2 


0.77 


1.2 


1.1 


0.86 


1.1 


1.1 


1.1 




d 


1.9 


0.56 


0.52 


1.5 


0.67 


0.66 


1.0 


1.0 


1.0 




e 


1.3 


0.78 


0.78 


1.2 


0.86 


0.86 


0.99 


0.99 


0.99 




f 


2.3 


0.53 


0.44 


1.7 


0.61 


0.56 


1.1 


1.1 


1.1 



Table B2A: Ratios of the rate of proton decay that would occur if chargino diagrams 
contributed only versus total proton decay rate for the three most dominant decay 
modes, for various values of the GUT scale parameters, when the ds operator is not 
included. 
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run 


m 










pchargino 






pcliargino 




no. 


d 
1 




ptotal 






ptotal 






ptotal 








max 


j3 = —a 


min 


max 


P = -a 


min 


max 


l3 = —a 


min 




a 


1.2 


1.1 


0.85 


1.1 


1.1 


0.87 


1.2 


1.1 


1.0 




b 


1.2 


1.2 


0.86 


1.2 


1.2 


0.86 


1.1 


1.1 


1.1 


I 


c 


1.1 


1.1 


0.87 


1.1 


1.1 


0.89 


1.2 


0.95 


0.92 




d 


1.2 


0.92 


0.92 


1.2 


0.92 


0.91 


1.1 


1.1 


1.1 




e 


1.1 


0.91 


0.91 


1.1 


0.91 


0.91 


1.2 


1.2 


1.1 




f 


1.2 


0.98 


0.96 


1.2 


0.97 


0.94 


1.1 


1.1 


1.1 




a 


1.4 


1.3 


0.70 


1.4 


1.3 


0.72 


1.4 


1.1 


1.1 




b 


1.4 


1.4 


0.76 


1.4 


1.4 


0.75 


1.3 


1.1 


1.1 


II 


c 


1.3 


1.3 


0.71 


1.3 


1.2 


0.75 


1.4 


0.92 


0.88 




d 


1.4 


0.89 


0.88 


1.4 


0.88 


0.86 


1.3 


1.3 


1.2 




e 


1.3 


0.81 


0.81 


1.3 


0.81 


0.81 


1.4 


1.4 


1.2 




f 


1.4 


0.99 


0.97 


1.4 


0.97 


0.93 


1.3 


1.3 


1.2 




a 


1.6 


1.4 


0.64 


1.5 


1.3 


0.69 


3.1 


0.53 


0.43 




b 


2.1 


2.1 


0.43 


2.0 


2.0 


0.47 


1.9 


0.64 


0.63 


iii(i) 


c 


1.4 


1.3 


0.72 


1.3 


1.2 


0.76 


3.1 


0.35 


0.31 




d 


2.2 


0.49 


0.47 


2.0 


0.53 


0.50 


2.1 


2.1 


0.87 




e 


1.4 


0.75 


0.75 


1.3 


0.76 


0.76 


3.1 


2.9 


0.65 




f 


2.5 


0.46 


0.41 


2.3 


0.50 


0.42 


1.8 


1.8 


0.94 




a 


1.3 


1.3 


0.74 


1.3 


1.2 


0.78 


3.1 


0.48 


0.38 




b 


1.7 


1.7 


0.55 


1.6 


1.6 


0.59 


2.2 


0.52 


0.51 


111(2) 


c 


1.2 


1.2 


0.81 


1.2 


1.2 


0.83 


2.8 


0.40 


0.36 




d 


1.7 


0.58 


0.57 


1.6 


0.62 


0.60 


2.4 


2.4 


0.72 




e 


1.2 


0.83 


0.82 


1.2 


0.83 


0.83 


2.7 


2.6 


0.55 




f 


2.1 


0.52 


0.47 


1.9 


0.56 


0.50 


2.1 


2.0 


0.82 




a 


1.3 


1.2 


0.78 


1.2 


1.2 


0.81 


2.9 


0.49 


0.39 




b 


1.6 


1.6 


0.61 


1.5 


1.5 


0.65 


2.4 


0.47 


0.46 


111(3) 


c 


1.2 


1.2 


0.84 


1.2 


1.1 


0.86 


2.4 


0.45 


0.41 




d 


1.6 


0.64 


0.63 


1.5 


0.67 


0.66 


2.5 


2.5 


0.64 




e 


1.2 


0.86 


0.85 


1.2 


0.86 


0.86 


2.4 


2.3 


0.54 




f 


1.9 


0.57 


0.53 


1.7 


0.61 


0.56 


2.2 


2.1 


0.75 



Table B2B: Ratios of the rate of neutron decay that would occur if chargino diagrams 
contributed only versus total neutron decay rate for various values of the GUT scale 
parameters, when the O13 operator is not included. 
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run 


m 


/pLLRR 
/ p^K + -, 








no. 


d 


\ r 


LLLL 




\ r 


LLLL 


















Vr 




a 


U.UUUU44 


U.Ull 


4.2 


U.UUUU / 6 


U.Ulo 


8.4 




b 


U.UUUUzz 


n nn/io 
U.UU4y 


1.9 


U.UUUUo 1 


U.UUoo 


3.6 


I 


c 


U.UUUUdo 


U.Ul / 


5.7 


U.UUUii 


U.Uzy 


12. 




d 


U.UUUUzz 


U.UU4o 


1.7 


U.UUUUo / 


U.UU / 1 


2.8 




e 


n nnnnR9 

U.UUUUOZ 


n nil 

U.Ull 


4.2 


n nnm n 

U.UUUIU 


U.Ulo 


6.5 




f 


U.UUUUid 


U.UUol 


1.3 


U.UUUUzd 


u.uuou 


2.1 




a 


U.UUUU4(5 


U.Ull 


5.1 


U.UUUU <o 


U.Ulo 


10. 




b 


U.UUUUzo 


U.UU4y 


2.6 


U.UUUUoo 


U.UUoo 


5.4 


II 


c 


U.UUUU / i 


U.Ul / 


5.7 


U.UUUiz 


n noo 
U.Uzy 


11. 




d 


U.UUUUzz 


U.UU44 


2.2 


U.UUUUo / 


U.UU / o 


3.5 




e 


nofififii 

U.UUUUUl 


U.Ull 


4.2 


n nnni n 

U.UUUIU 


U 01 Q 
u.uiy 


6.4 




f 


U.UUUUiO 


U.UUol 


1.7 


U.UUUUzO 


U.UUol 


2.7 




a 


U.UUUiy 


U.Uoz 


6.4 


U.UUUoi 


U.Uoo 


10. 




b 


U.UUUU 1 i 


U.Ulo 


3.4 


U.UUUiz 


U.Uzo 


5.8 


iii(i) 


c 


U.UUUoU 


U.U4o 


7.6 


U.UUU4y 


U.U {\) 


12. 




d 


u.uuuuoy 


U.U14 


3.8 


u.uuuuyo 


v.vZZ 


6.8 




e 




049 


8.5 


n unn9^i 


Ofi4 


16. 




f 


0.000042 


0.0093 


2.8 


0.000069 


0.015 


4.8 




a 


0.00029 


0.048 


10. 


0.00046 


0.081 


16. 




b 


0.00011 


0.020 


5.3 


0.00018 


0.034 


9.1 


111(2) 


c 


0.00044 


0.066 


12. 


0.00072 


0.12 


19. 




d 


0.000088 


0.020 


6.1 


0.00014 


0.032 


11. 




e 


0.00023 


0.061 


13. 


0.00037 


0.095 


24. 




f 


0.000062 


0.014 


4.4 


0.00010 


0.022 


7.6 




a 


0.00035 


0.058 


13. 


0.00056 


0.099 


20. 




b 


0.00013 


0.024 


6.8 


0.00022 


0.042 


12. 


111(3) 


c 


0.00054 


0.081 


15. 


0.00088 


0.14 


24. 




d 


0.00011 


0.024 


7.8 


0.00018 


0.040 


14. 




e 


0.00028 


0.075 


16. 


0.00046 


0.12 


30. 




f 


0.000076 


0.017 


5.6 


0.00012 


0.027 


9.7 



Table B4: Ratios of the rate of proton decay that would occur if LLLL operators 
contributed only versus the rate of proton decay that would occur if LLRR operators 
contributed only, for each of the three anti-neutrino generations, for various values of 
the GUT scale parameters, when the C13 operator is not included. 
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run 
no. 



\ 



LLLL 

LLHH 



LLHk 



/f+- 



10-6 
10-*^ 
10-6 

10^6 
10-6 



P —• 

K + -U 



\ 



LLltJt 



a 
6 

c 
d 

e 

f 



0.062 
0.13 
0.043 
0.12 
0.045 
0.17 



2.8 X 
2.8 X 
2.8 X 
2.8 X 
2.8 X 
2.7 X 



0.78 
1.7 

0.49 
1.9 

0.80 
2.7 



10-6 
10-6 
10-6 
10-6 
10-6 



0.0084 
0.0085 
0.0085 
0.0084 
0.0085 
0.0082 



0.24 
0.53 
0.18 
0.57 
0.24 
0.79 



II 



a 
b 
c 
d 
e 
f 



0.071 
0.15 
0.049 
0.15 
0.057 
0.22 



3.4 X 

3.5 X 
3.5 X 

3.4 X 

3.5 X 
3.3 X 



0.95 
2.1 
0.61 
2.3 
0.94 
3.2 



10^ 
10-6 

10-6 
10-6 
10-6 
10-6 



0.010 
0.011 
0.010 
0.010 
0.011 
0.010 



0.20 
0.38 
0.18 
0.45 
0.24 
0.60 



III(l) 



a 
b 
c 
d 
e 
f 



0.015 
0.043 
0.010 
0.050 
0.020 
0.069 



3.0 X 
3.0 X 
3.0 X 
3.0 X 
3.0 X 
2.9 X 



0.28 
0.69 
0.21 
0.67 
0.22 
0.95 



10-6 

10-6 
10-6 
10-6 
10-6 



0.0091 
0.0092 
0.0092 
0.0091 
0.0092 
0.0089 



0.16 
0.30 
0.13 
0.26 
0.12 
0.36 



111(2) 



a 
b 
c 
d 

e 
f 



0.010 

0.029 

0.0068 

0.034 

0.013 

0.047 



3.0 X 
3.0 X 
3.0 X 
3.0 X 
3.0 X 
2.9 X 



0.19 
0.46 
0.14 
0.45 
0.15 
0.64 



10-6 

10-6 
10-6 
10-6 
10-6 



0.0090 
0.0091 
0.0091 
0.0090 
0.0091 
0.0088 



0.10 
0.19 
0.085 
0.16 
0.076 
0.23 



111(3) 



a 
b 
c 
d 
e 
f 



0.0086 

0.023 

0.0056 

0.028 

0.011 

0.038 



3.0 X 
3.0 X 
3.0 X 
3.0 X 
3.0 X 
2.9 X 



0.16 
0.38 
0.11 
0.37 
0.12 
0.53 



0.0091 
0.0092 
0.0091 
0.0090 
0.0092 
0.0089 



0.080 
0.15 
0.068 
0.13 
0.061 
0.18 



Table B5A: Ratios of partial decay rates for p — > K'^V, which compare the importance 
of the LLLL and LLRR operators for each generation of anti-ucutrino versus contri- 
bution of the LLRR operator of the third generation anti-neutrino for various values 
of the GUT scale parameters, when the ds operator is not included. 



23 



run 
no. 



II 



iii(i) 



111(2) 



111(3) 



a 
b 

c 
d 
e 
f 



a 
b 
c 
d 
e 
/ 



a 
b 
c 
d 

e 
f 



a 
b 
c 
d 

e 
f 



a 
b 
c 
d 
e 

1 



\ 



LLLL 
7r + -e 

LLliH 



LLRR 

LLHH 



\ 



LLUll 



0.038 
0.076 
0.026 
0.076 
0.027 
0.11 



2.8 X 
2.8 X 
2.8 X 
2.8 X 
2.8 X 
2.7 X 



10"'' 
lO"'^ 
10"*^ 

10-6 

10^6 
10-*^ 



0.48 
1.0 

0.30 
1.2 

0.48 
1.7 



0.0085 
0.0085 
0.0085 
0.0085 
0.0085 
0.0083 



0.044 
0.092 
0.030 
0.095 
0.034 
0.13 



3.5 X 
3.5 X 
3.5 X 
3.5 X 
3.5 X 
3.4 X 



10^ 
10"*^ 
10-6 
10-6 

10-6 
10-6 



0.58 
1.3 

0.37 
1.4 

0.57 
2.0 



0.010 
0.011 
0.011 
0.010 
0.011 
0.010 



0.0097 
0.025 
0.0063 
0.031 
0.012 
0.043 



3.0 X 
3.0 X 
3.0 X 

3.0 X 

3.1 X 
3.0 X 



10^ 

10-6 
10-6 
10-6 
10-6 
10-6 



0.17 
0.40 
0.12 
0.42 
0.14 
0.60 



0.0092 
0.0092 
0.0092 
0.0092 
0.0093 
0.0090 



0.0065 
0.017 
0.0042 
0.021 
0.0081 
0.029 



3.0 X 
3.0 X 
3.0 X 
3.0 X 
3.0 X 
2.9 X 



10^ 

10-6 
10-6 
10-6 
10-6 
10-6 



0.11 

0.27 

0.078 

0.28 

0.096 

0.40 



0.0091 
0.0091 
0.0091 
0.0091 
0.0092 
0.0089 



0.0054 
0.014 
0.0034 
0.017 
0.0066 
0.024 



3.0 X 
3.0 X 
3.0 X 
3.0 X 
3.0 X 
3.0 X 



10-6 
10-6 
10-6 
10-6 
10-6 
10-6 



0.092 
0.22 
0.064 
0.23 
0.079 
0.33 



0.0091 
0.0092 
0.0092 
0.0091 
0.0092 
0.0090 



LLLL 

LLltJt 



0.12 

0.28 

0.084 

0.35 

0.15 

0.48 



0.095 
0.19 
0.088 
0.29 
0.16 
0.38 



0.096 
0.17 
0.080 
0.15 
0.064 
0.21 



0.062 

0.11 

0.052 

0.093 

0.042 

0.13 



0.049 
0.085 
0.042 
0.072 
0.034 
0.10 



Table B5B: Ratios of partial decay rates for p — > tt+I/, which compare the importance of 
the LLLL and LLRR operators for each generation of anti-ncutrino versus contribution 
of the LLRR operator of the third generation anti-neutrino for various values of the 
GUT scale parameters, when the O13 operator is not included. 



24 



run 
no. 


I 


II 


III(l) 


111(2) 


111(3) 




24.43 


24.36 


24.51 


24.65 


24.75 


Mqut 


2.498 X 10^^ 


3.172 X 10^*^ 


3.327 X 10^^ 


2.857 X 10^^ 


2.513 X 10^^ 


£3 


-0.04760 


-0.04886 


-0.04342 


-0.04420 


-0.04550 


A 


0.7640 


0.8067 


0.8523 


0.8867 


0.8872 


B 


0.05798 


0.06019 


0.06254 


0.06533 


0.06607 


C 


0.00008824 


0.00009204 


0.00009550 


0.00009801 


0.00009809 


E 


0.01063 


0.01111 


0.01154 


0.01182 


0.01180 




1.762 


1.765 


1.767 


1.765 


1.763 


tan (3 


52.71 


54.31 


55.32 


55.79 


55.87 




80.0 


80.0 


160. 


240. 


300. 




280. 


240. 


170. 


170. 


170. 


Too 


400. 


700. 


1400. 


1400. 


1400. 


mn^ 


706.3 


994.4 


1858. 


1859. 


1855. 




635.9 


865.6 


1599. 


1592. 


1585. 


Ao 


322.2 


458.4 


-982.4 


-1079. 


-1274. 



Table B6: Values of the GUT scale parameters used in Tables Bl through B5. All 
dimensions in GeV units. 



5 Discussion of the results 

5.1 Overall Rates 

Many significant results appear from the tables of the previous section. First, 
comparing the rates of proton decay predicted in the tables above with the 
results of experimental searches for proton and baryon- number violating neutron 
decay summarized in Table ^, it can be seen that the predicted upper bounds 
on the lifetimes for nucleon decay are above, and, in most cases, well above, the 
experimental lower bounds. The loop integral /(M„atioino, , M sg^ark ) 

goes roughly like 1/Af^,,„„^fc in the limit where squarks and sleptons are much 

heavier than gauginos. Hence, the decay rates go naively like {m^^^ + /i^)/TOo, 
where iir = ii{Mz). This approximation roughly explains the dependence of 
the nucleon decay rates on toq, ii{Mz), and mi/2 seen in Tables AlA, AlB, 
BIA, and BIB. 

5.2 LLRR vs. LLLL Operators 

Secondly, LLRR operators dominate over LLLL operators for the third genera- 
tion anti-neutrino, for the decays into KV and ttV. (See Tables A4 and B4.) We 
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Table 4: Current experimental lower bounds on the various partial lifetimes of 
the nucleons 



rip - 


K+V) 


> 


1.0 X 10^2 


yrs 


t{p - 




> 


.25 X 10^2 


yrs 


t{p 




> 


2.7 X 10^2 


yrs 


t{p - 




> 


.69 X 10^2 


yrs 


rip 




> 


5.5 X 10^2 


yrs 


t{p - 


rye^) 


> 


1.4 X 10^2 


yrs 


T{n — 




> 


.86 X 10^2 


yrs 


T(n - 




> 


1.0 X 10^2 


yrs 


T{n 


rfV) 


> 


.54 X 1032 


yrs 


T(n 




> 


1.0 X 10^2 


yrs 


T{n — > 




> 


1.3 X 10^2 


yrs 



can gain an intuitive understanding of why this is if we neglect the gluino contri- 
bution to the decay rates and look at approximate formulas for the rate of decay 
due to charginos. The loop integral I{a, b, c) is a relatively smooth function of 
the masses and as a result, to a very good approximation, when calculating 
chargino diagrams for the third generation anti-neutrino, sums over gamma 
matrices of the form Tj^ \jJ^*R xj^i^x^ ^7 c) are approximately zero while sums of 
the form aj^l Aj"-^(^A: b, c) and F^j xJ^r \j^i^>^^ ^' ^'^^ approximately equal 
to 6ijI{Vti^ , 6, c) and Sijl{fli^, b, c), respectively, due to the orthogonality of the 
gamma matrices. Hence, 

Wn^Mt i^U,RYuVKM)xkTlj^R Xk'(XE,RYe)p3^*E,Rpl'Cla^^^cl^'' U+ 2„J7_ 2«TOx./(x«, "A, Cp) 

~ 16^2M -^r At ( M ) 3fc C* ^ ^^^U+ 2nU - 2n»™x™ ^(Xn , , 63^ ) 

_ ' . . . . . 

See fig. y for the Feynman diagram giving the dominant contribution to eqn. 



167r2M ln"Zx„{52(^ifM)»'i(V"^M)i'lCgg'''^'cg;''l^/(Xn,X^ 
+5|(V^ifM)fe'fcC,qf^^C,z'='l3/(x„,Mfe^,g3i)} (4) 

The loop integral factors 'Ylin^+ o-nU- a'TL^ixnliXnTb, c) can further be ap- 
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Figure 2: Feynman diagram that gives the dominant contribution to C{-^^ 



{ud){diy) 



> 



\ 




proximated 

VC/^ U : m- I{y be) ^ \ My,inoI{W±,b,c) if a = a' = 1 



fijiI{H±,b, c) ii a = a' = 2' 

Hence, 

(5) 

The integral I{a,b,c) is approximately log(6^/c^)/(6^ — c^) in the limit where 
a b,c, i.e. when squarks and sleptons are much more massive than gauginos, 
which generally is the limit we are interested in. Using the fact that the first and 
second generation squarks are approximately degenerate, and that the first and 
second generation squarks are usually more massive than the third generation 
squarks and sleptons, the ratio can be further approximated 

where m^'^^^^i^ is the mass of the first and second generation squarks. 

Several factors contribute to the fact that this ratio is often greater than 
one. Since the third generation squarks and sleptons are lighter than the first 
and second generation squarks, the ratio is significantly enhanced by the factor 
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K'quark)'log("^l3„/"^i3«)/("^i3^. -"^kj ~ K'quark)V(max(mi3^,mg3^))2. 
Second, the ratio 

is itself typically of order unity for many values of the GUT scale initial parame- 
ters. Note, this ratio is greatly enhanced in the regime of large tan /3 considered 
in this paper. 

Finally, the ratio ^irI Mwino plays a critical role in whether the LLRR op- 
erators dominate over LLLL operators in the third generation and whether the 
third generation anti-neutrino dominates over the second-generation. Compar- 
ing Tables A4, ASA, A5B, B4, B5A, and B5B with Tables A6 and B6, we see 
that there is a direct correlation between /ifi/mi/2 and the dominance of the 
third generation LLRR operators. When /Ltij/mi/2 is small, the third genera- 
tion LLRR operators can be suppressed. Moreover, because the third generation 
LLRR operators are suppressed when /ij?/mi/2 is small, the second generation 
anti-neutrino contributes more significantly than the third when is 
small.|^ In particular, in runs I and II of the A and B tables, /J.i?/TOi/2 is small 
(~ .3) while in runs III(l), 111(2), and 111(3), /i_R/?7ii/2 is near 1 or greater, and 
Mi?/™i/2 increases as one goes from run III(l) to run 111(2) to run 111(3). As 
a result, looking at the tau anti-neutrino columns of Tables A4 and B4, we see 
that, for any particular model, the entries of those columns for runs I and II 
are smaller than they are for runs III(l), 111(2), and 111(3), and that the entries 
in those columns increase steadily in going from run III(l) to 111(2) to 111(3). 
Similarly, looking at the 3rd columns of Tables A5A, A5B, B5A, and B5B, we 
see that the LLLL operators of the second generation anti-neutrino are fairly 
significant to the overall decay rate in runs I and II, while they are not quite 
as significant in runs III(l), 111(2), and 111(3), and that the significance of the 
LLLL operators of the second generation anti-neutrino continually decreases in 
going from run III(l) to 111(2) to 111(3). 

For the first and second generation anti-neutrinos, on the other hand, the 
LLRR operators are negligible because they are suppressed in comparison to the 
LLLL operators by the up and charm Yukawa couplings, respectively. Thus, the 
entries in the electron and muon anti-neutrino columns of Tables A4 and B4 are 
fairly small, and the entries in columns 2 and 4 of Tables A5A, A5B, B5A, and 
B5B are fairly small. In comparison, the second generation anti-neutrino LLLL 
operator is the most significant of the LLLL operators, but the third generation 
LLLL operator is not negligible in comparison to the second generation LLLL 
operator. (See columns 1, 3, and 5 of Tables ASA, A5B, B5A, and B5B.). 



''Arnowitt, ct al. observed that LLRR operators can be significant to nucleon decay rates 
under certain circumstances in ref. |20l . 
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5.3 



p — > TT^Z/ VS. p 



Secondly, we see that under certain circumstances, the decay p n^V domi- 
nates over K'^V when the ©13 operator is included in model 4(c). When the 
third generation anti-neutrino dominates over the the other generations, we have 
the approximate result 



T{p^ K+v) 



3.9- 



(Vft:M)3ic* 



[llg* 3]3 



3.9 



iVKM)31 



KM }32 



3] 3 



[12 



{Vkm)32 ud ue 



C* 3]3 _|_ ^* [llg* 3]3 



(8) 



In Appendix 4, we show that \{Vkm)3i\ = \Vtd\ increases when the O13 
operator is included. The increase in the ratio of the rate of p ^ tt^V versus p 
K^v when the ©13 operator is included can be attributed in large part to this 
increase in \Vtd\, since eqn. ^ contains a multiplicative factor of |Vtd/Vtsp. This 
increase is further enhanced by the fact that {Vkm)3ic^J'^^(Z.J^^^ has roughly 

the opposite sign of (l^ifA£)32C*^'^^c*g^l^, and hence increasing \Vtd\ decreases 
the denominator in eqn. Q. Thus, the addition of the O13 operator in model 
4(c) increases the ratio of the p tt+F decay rate to p ^ K^V, provided that 
the third generation anti-neutrino dominates. 

Whether the third generation dominates over the second depends on the 
ratio of Thus, p tt^V will be larger than p —>■ K^V if the O13 

operator is included and /ii?/TOi/2 is not much smaller than one. Thus, in runs 
III(l), 111(2), and 111(3) of Table Al, iJLR/mi/2 is approximately one or bigger, 
and as a result the third generation anti-neutrino dominates the rate of decay, 
and the ratio of the rate of decay into tt+F versus the rate of decay into K'^V 
is significantly enhanced in comparison to the runs without the ©13 operator 
in Table Bl. On the other hand, in runs I and II, Atij/TOi/2 is small, and as a 
result, the second generation dominates and the ratio of the rate of decay into 
TT+T? versus K^v remains near what it was without the O13 operator. 



5.4 "Generic" SU(5) vs Large tan/3 SO(IO) models 



5.4.1 71 — > ■k'^v vs. n 



Furthermore, the tables of the previous section show some important differ- 
ences between the nucleon decay predictions for our SO (10) model versus the 
predictions of a generic SUSY minimal SU(5) model. Because the effective color 
triplet Higgs mass is constrained to be lower than around 10^^ GeV in SUSY 
minimal SU(5) [Q, and because the lifetimes of the nucleons are proportional 



29 



to sin^ 2/3 m 



minimal SU(5) models use small tan f3 to be consistent 



with the experimental limits on proton decay. When tan/3 is small, LLRR op- 
erators can often be neglected. When LLRR operators are negligible, the ratio 
T(n —>■ TT'^17)/T{n — > rfp) just depends on chiral Lagrangian factors. 



T{n 7r"i7) 



2.8- 



J2- /3C("'')('^''*) + aC^"'')^''''*) 



P(j(ud){diy,) _ .i40aC("'')('*''') 



2.8 



(9) 



In contrast, when tan/3 is large, LLRR operators are not negligible. The contri- 
bution of LLRR operators to n rfU is significantly suppressed in comparison 
to its contribution to n ^ tt'^F by chiral Lagrangian factors. Hence, looking 
at Tables AlB and BIB, the rate oi n ^ tt^V can be anywhere from 2.9 to 
over 100 times larger than the rate of rt ^ 77F, depending on whether the third 
generation LLRR operators or the second generation LLLL operators dominate. 



5.4.2 



K v vs. p K^v 



Secondly, the ratio V{n K'^v)/T{p 
generic SU(5) models. Numerically, 



K^v) differs significantly from the 



T{p^ K+V) 



/3(1.14C("^)('^'^-) + 1.58C("'^)('*''')) +a(-.86C("")(''''') + 1.58C('"^)("''')) 



/3(.44C(«s)(di.,) _^ i.58C("'')(^''')) + a(.44C(«^)(''''») + 1.58C(«'^)('' 

(10) 

In the minimal SU(5) model, the c'^^'^^id-^i) operator is approximately equal 
to the C("'') m, |0[ |1, and, since LLRR operators are often negligible, 
r(n ^ K^V)/T{p K+V) K. 1.8. g]. However, in our SO(IO) model, LLRR 
operators tend to dominate. Not only are the chiral Lagrangian factors different 
when the LLRR operators dominate, but C''^"'*^^'*'''-' tends to point in the opposite 
direction as C'( Hence, Tin K°T')/T{p K+"U) is much larger than 
its value in the minimal SU(5) models. Indeed, n K^V can be over 18 times 
bigger than p — > K+v. 

Also noteworthy is the fact that when LLRR operators dominate, T{n ^ 
K'^V)/T(p K+V) is significantly higher when the Ois operator is included in 
comparison to when it is not. (For example, in run 111(3) V{n ^ K^v)/T{p 
K+V) is 18.1 with the ©13 operator included while it is no greater than 4.5 
for run 111(3) without the O13 operator.) Much of the enhancement can be 
explained by the fact that \Vtd\ is larger when the da operator is included. 
When LLRR operators dominate. 



.)) 



T{p^ K+V) 



— .86(7^"'*-'^'''''^ + 1.58C^"'^'('*''*^ 



(11) 
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Plugging cqn. ^ into this formula, this becomes 



r(n ^ K"V) 
T{p-> K+V) 



(12) 



Since iVKM)3icla^'^'^c*J^^ tends to have the opposite sign as (VKM)32c;rf'"c*/l3, 
T{n K'^v)/T{p K^v) is enhanced when \Vtd\ is increased. 



5.5 Sensitivity to "22" Clebsch 

Furthermore, by looking at Tables Bl through B6, we can determine how sensi- 
tive nucleon decay rate predictions are on the Clebsches that enter into the Cgq, 
Cqu Cud, and Cue matriccs. For each of the five runs of Tables Bl through B6, 
the only difference between the versions a through / in each run are the Uqq, 
UqU Hud and yue Clebsches of Table 1 that enter into the Cqq, Cqi, Cud, and Cue 
matrices. We see that the overall rate of decay is quite sensitive to the different 
choices for Clebsches. For example, with f3 — —a, t{p — > K^V) for run 1(e) is 
7 times larger than T{p K^V) for run 1(f). The branching ratios generally 
exhibit less sensitivity: V{n ^ 7r'^77)/r(n K'^V) exhibits virtually no sensitiv- 
ity and r (p — > tt+F) /F (p — > K^V) exhibits relatively mild sensitivity. However, 
branching ratios into less dominant decay modes can at certain times exhibit 
high sensitivity. For example, with /3 = —a, r(7i r]h')/T{n — > K^V) is over 12 
times larger for run 1(f) than it is for run 1(c). 



5.6 Gluino vs. Chargino contributions 

It can also be seen that the contributions of gluinos to the rate of nucleon 
decay is often not negligible. Indeed, in several examples, excluding the gluinos' 
contribution can lead to a decrease in the predicted rate of decay of greater than 
60%, in cases where gluinos constructively interfere, or an increase in the rate 
of decay by over 150%, where gluinos destructively interfere.^ Note also that 
whether gluinos constructively or destructively interfere depends heavily on the 
phase of the chiral Lagrangian parameter a,Tg{P/a). 



5.7 Proton decay from gauge boson exchange 

Finally we note that our analysis only includes the contribution to nucleon decay 
from the effective dimension 5 operators resulting from colored triplet Higgs 
exchanges. We have neglected the contribution to nucleon decay via heavy 
gauge boson exchange (effective dimension 6 operators) . This approximation is 
justified in our models for the dominant decay modes. For example, in order to 

*Goto, et al. observed that gluino loops can be important in the minimal SU(5) model in 
rcf. g]. 
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obtain the £3 of run 111(1), we can choose the vevs ai, 02, a, and a singlet field 
<S4, defined in paper I, which enter into the SO(IO) breaking sector of the theory, 
to be 2.0 X 10^6 GeV, 1.0 x lO^^ GeV, 6.0 x lO^^ GeV, and .66 x 10^^ GeV, 
respectively, i.e. all of order Mqut- Then the masses of the gauge bosons 
contained in 5'O(10)/(S't/(3) x SU{2) x U{1)) - X±, Q±, and are 
3.2 X 10^6 GeV, 1.2 x 10^^ GeV, 1.3 x lO" GeV, and 1.3 x lO" GeV, respectivelyg. 
The decay mode which would be the most dominant if all other contributions 
to proton decay except the contribution due to gauge boson exchanges were 
neglected is p ^ Tr^e"*". With the above gauge boson masses, the partial proton 
lifetime due to heavy gauge boson exchanges for p ir'^e'^ is 1.2 x 10'^* yrs, 
corresponding to a branching ratio of order < 10^**. Of the decay modes hsted 
in Tables AlA and BlA, gauge exchange is competitive only with p r/fj.^. 



Conclusions 

We have shown in this paper that model 4(c) of paper I predicts nucleon decay 
rates consistent with all current experimental bounds, while using values of 
GUT parameters that give fermion masses, mixing angles, and gauge couplings 
in good agreement with experimental observations. Our main results can be 
found in Tables AlA, AlB, and A6. We conclude that our model predicts that 
nucleon decay is likely to be observed by SuperKAMIOKANDE or ICARUS, 
which are expected to probe nucleon lifetimes up to around 10'^'* yrs j2^, p^ , 
for various decay modes predicted by GUTs. 

In order to avoid this conclusion one would need to make squarks and slep- 
tons "unnaturally" heavy and beyond the reach of LHC or increase the effective 
color triplet Higgs mass Aft, which would require a supermassive Higgs doublet 
in the GUT desert with mass many orders of magnitude lower than the GUT 
scale and at least an order of magnitude lighter than any other particle getting 
mass around the GUT scale. Moreover, we have chosen the poorly known chi- 
ral Lagrangian parameter to be at the lowest value suggested by the data. 
If it could be shown that \(3\ lies in the higher range of its current bounds, 
non-observation of nucleon decay by SuperKAMIOKANDE and ICARUS could 
make our model unnatural for any reasonable values of the squark and slepton 
masses.]^ For these reasons, we believe that these models, if correct, necessarily 
lead to observable nucleon decay rates. 

We have shown that LLRR operators are not only significant, but often 
dominate, nucleon decay in the large tan /3 regime - as long as fJ-R/mi/2 is 

^Note, the gauge boson is the massive gauge boson from the 24 representation of 
SU(5); is the gauge boson from the 10, 15, 10, and 15 representations of SU(5) which 
is in the (3,2, i) and (3, 2, -i) representations of SU(3)xSU(2)xU(l); and [/± and £;± are 
gauge bosons from the 10 and 10 representations of SU(5) which are in {(3, 1, (3, 1, — ^)} 
and {(1, 1, -2), (1, 1,2)} representations of SU(3)xSU(2)xU(l), respectively. 

Recall, however, that the chiral Lagrangian approach tends to overestimate nucleon decay 
rates llql and thus underestimates the lifetimes. 
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not very small. As a result, if nucleon decay is observed, there are two key 
experimental observables that may distinguish between large or small tan /? 
SUSY GUTs: the ratios r(n ^ K^17)/T{p K+v) and T(n ^ riv)/T{n 
tt^F). In particular, evidence for a neutron lifetime (| — ^) x the proton lifetime 
would be a strong indication for large tan (5 SUSY GUTs. Observation of 
Tin ^ rjU) /r [n tt'^V) significantly lower than than the predicted value when 
LLRR operators are negligible could also indicate large tan (3 SUSY GUTs. 

We have also shown that gaugino loops cannot be neglected when calculating 
proton decay rates in models such as ours. In fact, neglecting gaugino loops 
could lead to an underestimation of the decay rates of over 60% or overestimation 
of over 150%. 

Finally, we have studied the sensitivity of nucleon lifetime and branching 
ratio predictions on the "quality" of the predictions that these models make for 
fermion masses and mixing angles. In the models we have analyzed, the entries 
in the Cqq, Cqi, Cud, and c„e matrices are related to entries in the y„, Yd, and Yp 
matrices by Clebsches which depend on the version of the model being consid- 
ered. As we have seen, the lifetimes and branching ratios can be quite sensitive 
to the choice of Clebsches — some predictions vary by nearly an order of magni- 
tude depending on the choice of Clebsches. In models 4(a) through (f), without 
the Oi3 operator, the different Clebsches have no effect on the predictions for 
fermion masses and mixing angles. Comparing models 4(a) through (f), without 
the 013 operator, which is consistent with fermion masses and mixing angles at 
2(7, with model 4(c), with the C13 operator, which is consistent within la, one 
is lead to conclude that fitting the data within 1 or 2 cr can have a significant 
effect on the nucleon decay predictions. Thus "predictions" for nucleon decay 
lifetimes and branching ratios cannot be expected to be any better than the 
complementary predictions for fermion masses and mixing angles. 
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Appendix 1: How the Dimopoulos-Wilczek mech- 
anism can produce Mt ^ Mqut 

Using the Dimopoulos-Wilczek mechanism for doublet-triplet splitting, the part 
of the superspace potential in our model giving doublet-triplet splitting is 

Wd-t = lOi^ilOa + lOi^*, 
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where Ai is a 45 representation getting a vev of order Mqut in the baryon 
minus lepton number direction, lOi and IO2 are 10 representations, and <S* is a 
singlet. The Higgs doublet and triplet mass matrices for our model are 



lOi 
IO2 



lOi 


-ai 



IO2 

ai 



(13) 



lOi IO2 

101 / 

102 I 5* 



(14) 



Note, the lOi field is the only 10 representation that couples to ordinary (Stan- 
dard Model) fcrmions. Thus, the lOi contains the two Higgses of the Mini- 
mal SUSY Standard Model [MSSM], and the baryon- number violating effec- 
tive operators obtained by integrating out the color triplets are proportional to 
1/ Mt = {M^^)ii = S^/af, the lOi, lOi entry of the inverse of the color triplet 
mass matrix. An effective color triplet mass of around 10^^ GeV is obtained 
with oi around the GUT scale, and 5* around 10^^ GeV. Thus, Alt ^ Mqjjt 
means that there is an electroweak doublet several orders of magnitude lighter 
than the GUT scale, not that there is an actual color triplet with mass greater 
than M Planck fH] . 

Appendix 2: Feynman diagrams 



U+ 2n)VKM]xj5'fj 



"■p 



X- 
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m^,, ijU+ln 



— »— 



U-2n] 



x+ 
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Appendix 3: Formulas for nucleon decay in terms 
of chiral Lagrangian factors 

Using the chiral Lagrangian techniques of ref. jl^, the rates of nucleon decay 
are the following. 



Tip ~> TT+i>i) = 
T{p 776+) = 



(mg - m^)^ 
327rm3/2 



3mB 



+ [pC^^d)(s^^) + aC('"^)("''')][l + -^{D + 3F)] 
^p(j(ud){d^,) ^ + D + F) 

1 



3{ml 



2 9\9 



» . ^il /3C("'')("^')[1 + -(3^^ - D)] - aC('"^)("^')-[l - (3F - D)] 
bATTf^m-^ 3 3 



^^(«d)(«e.)[i _^ 1(3^ _ - aC('"')("'=*)i[l - (3F - D)] 



+ , 

(mg - TO^)2 
327r/2m3 



} 



rriB 



TTlB 



niB 



+ [f3C^'^'>(^^ +aC^'''^^^^'>]{l + D + F)^} 



TUB 
1 



} 



2 ^2 



{ml - m|-) 



3m B 3m B 



+ {f3C<^^d)(s^.^) + aC^"'^) ))(1 + + 3F)) 

3m B 
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r(n^7r°I7,) 



" a2 

647r/2 ^ 



P(j{ud){dui) ij^(j{ud){dui) 

■i{ml - mlf 



{\ + D + Ff 



All 



p(j{ud){d^,) + -{3F - D j) - aC("'^)(''''')i(l + £» - 3F) 
3 3 



^^(u(i)(«ei) _|_ ^^(«d)(«ei) 



/3C 



Kl + I^ + F)^ 



where is an average Baryon mass satisfying ttib ~ m^: ~ ttla and all other 
notation follows Here, all coefficients of four-fermion operators are eval- 

uated at Mz- Al takes into account renormalization from Mz to 1 GeV, and 



is approximately equal to .22. 
grangian formulas given in ref. 



we take D = .81, F = .44 [|, and = 139 MeV || 



These formulas reduce to the chiral La- 
for /? 7^ when a ~ Q. In the calculations. 



Appendix 4: Why \Vtd\ increases when the O13 is 
included in model 4(c) 

It can be seen that our Yd and Ye matrices, which have the general form. 




zC 
x'B 




with C, D B, E A, can be diagonalized by multiplying the matrices on the 
left and right by matrices S and T, respectively, where 



( 



{zC- 



V 



1 

A ^ ' yE 



-{zC- 



A 

1 

xB 
A 



/ 1 



xuBDA5\ 
A ^ I 



A ' i/BA 

Our Yu matrix, which has the general form 



^ > yE 



xBe' 



yEA 



{zC- 



xu'BD „-iS\ e"^ 

'yE 



'-{zC- 

xB 
A 



X uBD A5 



1 

A 
x'B 
A 

1 
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{ud){dv)[G] 
ijkl 
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with C,D B,E A and C < xx'B'^/A, is diagonalized by the S and T 
matrices 



1 A I x'uBD „iS\ zC 

XX' A ^ ) ~FB 
a: I A I'^r^ x'uBD „-i5\ -1 

tt-D Q — iS xB 2 



1 ^(zC- ^lii^e-^'') li^e-*'^ 



Therefore, at Mgut 



A AE 2A3^A'e 

« -12e-n|g + i^e^(^+«) (15) 

Since S + (f) fa 30°, the De^^^'^'^^ /{'oA) term will increase | Vtc(| provided that the 
BC /{AE) term does not decrease too much when the Oia operator is included. 

What effect does the ©13 operator have on the A, B, C, and E parameters? 
The A parameter will not be affected at all because A is fixed purely by the 
third generation masses mi, and rrir 0. Moreover, the first generation Yukawa 
coupling Ag of the Ye matrix is approximately equal to 



243 



109 „„B 
27 - " A^^' 



(72 CD—"'^ 



(16) 



at Mgut- Since 6 is typically close to 2tt, the term of the absolute value linear 
in D decreases the eigenvalue. Therefore, C must be increased to compensate 
for the Oi3 operator. 

On the other hand, when B is evaluated by using the global analysis of 
ref. B actually decreases when the O13 operator is included. It will be shown 
in ref. that when a global analysis is done, Bk, the bag constant which 
comes into the theoretical formula for the experimental observable measuring 
CP violation, and \Vcb\ come out too high; and iKib/Vcfcj comes out too low, 
in comparison with their experimental measurements, for model 4 without an 
©13 operator. Since Vcb ~ C\^d ~ Xu\B/A where C is a renormalization group 
factor, this means that B is too high. It will also be shown in ref. |^ that 
unless an 013 operator is included, Vcb and Vub/Vcb cannot be corrected by 
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changing the parameters without further increasing Bk, which is ah'cady too 
high. Furthermore, B and E are related by the equation 

\3AEe"'' -B^\ « A^A^. 

Since Re (/> < 0, is lower when the O13 operator is not included because B is 
higher than it should be. 

However, when the O13 operator of paper I is added to model 4(c), it is 
possible to lower B to bring it in line with the experimental value for \Vcb\ 
while simultaneously having reasonable values for Bk as well as for the other 
observables [1). The net effect of including the ©13 operator to model 4(c) is 
that B is typically decreased by ~ 10%, E is increased by ~ 18%, and C is 
increased by ~ 25%. Therefore, BC/{AE) actually decreases slightly (~ 3%). 
However, the decrease in that term is more than compensated for by the increase 
in \Vtd\ due to the ^e**^''+"^^ term. The net result of the inclusion of the ©13 
operator is an increase in \Vtd\ typically of ~ 11%. 
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